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We investigate few body physics in a cold atomic system with synthetic dimensions (Celi et ah, 
PRL 112, 043001 (2014)) which realizes a Hofstadter model with long-ranged interactions along 
the synthetic dimension. We show that the problem can be mapped to a system of particles (with 
SU(M) symmetric interactions) which experience an SU(M) Zeeman held at each lattice site and 
a non-Abelian SU(M) gauge potential that affects their hopping from one site to another. This 
mapping brings out the possibility of generating non-local interactions (interaction between particles 
at different physical sites). It also shows that the non-Abelian gauge field, which induces a flavor- 
orbital coupling, mitigates the “baryon breaking” effects of the Zeeman field. For M particles, the 
SU(M) singlet baryon which is site localized, is “deformed” to be a nonlocal object (“squished” 
baryon) by the combination of the Zeeman and the non-Abelian gauge potential, an effect that we 
conclusively demonstrate by analytical arguments and exact (numerical) diagonalization studies. 
These results not only promise a rich phase diagram in the many body setting, but also suggests 
possibility of using cold atom systems to address problems that are inconceivable in traditional 
condensed matter systems. As an example, we show that the system can be adapted to realize 
Hamiltonians akin to the SU (M) random flux model. 

PACS numbers: 37.10.Jk, 21.45.-v, 67.85.-d 


Emulation of quantum systems of interest to a wide 
range of physics from condensed matter, high energy field 
theory etc., is made possible with cold atoms.pQ Recent 
experimental developments in the field stand testimony 
to this. nu Adding to the excitement and soaring in¬ 
terest in this area is the possibility of realizing systems 
that have novel physics and yet exceedingly difficult or 
not realizable in conventional condensed matter systems. 
Systems with SU ( M ) (M > 2) symmetries are one such 
example. Indeed, SU(M ) symmetric spin models have 
interesting phases and phase transitions, jSHU] as do 
Hubbard models with SU(M) symmetry. [TTH20] Several 
theoretical [2TH28] and experimental [291 - 135] works have 
explored physics of cold atomic systems with SU(M ) 
symmetry. Some of experimentally realized SU (M) sys¬ 
tems are 6 Li(M = 4) [36], 173 Yb ( M = 6)[32j [35j [37] . 
and S7 Sr(M = 10). 03 EH ES 

Celi et al m proposed the concept of “synthetic di¬ 
mensions” which achieves the goal of realizing finite sized 
“strip” of a Hofstadter model. Their idea, illustrated in 
Fig. 0 involves atoms with M internal states (labeled 
1... 7 ) in a ID (this can also be in higher dimensions) 
optical lattice. The hopping of the atoms from a site j 
(with coordinate Xj = jd, d is the spacing of the opti¬ 
cal lattice) to its neighbour does not change its internal 
state, and the amplitude t is independent of 7 . The inter¬ 
nal states at a site j are now coherently coupled such that 
an atom in state 7 at site j can “hop” to the state 7 + 1 
at j with an amplitude This produces, as shown 
in Fig. [lj a square lattice strip of finite width with M 
sites along the “synthetic dimension”. Since the coher¬ 
ent coupling is produced by a light of wavenumber kg, we 
have CP = Cl 1 e ~ lklXj , and this results in an atom picking 


up a phase factor e~ lked upon hopping around a plaque- 
tte. Choosing kgd = 27r| where p and q are relative 
prime integers, provides a realization of a finite strip of 
the Hofstadter model [39] with a p/q flux per plaquette. 
Very recent experimental realization[37] of this scheme 
bolsters the possibilities and scope of this research direc¬ 
tion. 

Another interesting aspect of the problem is that 
the SU(M) symmetric interactions between the atoms 
at a site j manifest as “infinite-ranged” (distance- 
independent) interactions along the synthetic dimension. 
For example, two atoms at site j (see Fig. [l]) with 7 = 1 
and 7 = 2 will interact with the same strength as 7 = 1 
and 7 = 4(M). It is the physics of such a system that 
is the subject of this paper, i. e., to understand inter¬ 
play between the flux p/q and the SU(M) interactions. 
It is essential to focus, as we do, on the physics of few 
particles since it provides crucial insights into construct¬ 
ing a many body phase diagram of the system. Pre¬ 
vious studies [23] i 25H27| of fermionic atoms with attrac¬ 
tive 5F(M)-interactions in a simple Id lattice (no flux, 
L e., | = 0, f2 7 = 0) show the existence of SU(M) sin¬ 
glet “baryons” and quasi-long-range color superfluidity of 
these baryons (see [40] for a review). The central question 
that we address in this paper is the fate of these baryons 
in the presence of flux in the synthetic dimension. 

Here we show that the Hofstadter model with infinite- 
ranged interaction along the synthetic dimension [38] can 
be mapped to a problem of M-flavor particles with 
SU ( M ) symmetric interactions hopping on the physical 
lattice with an on-site SU(M) Zeeman potential (de¬ 
termined by the couplings fl 7 ) along with a SU(M) 
gauge field (determined by the flux p/q , and f2 7 ) that 
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FIG. 1 . (Color online) Synthetic dimension: Atoms with 
M {Mm 4 is illustrated) internal states (labelled by 7 ) hop 
on a ID optical lattice (indicated by the alternating pink and 
cyan pattern) with amplitude t. The internal hyperhne states 
(indicated by thick brown lines with a red halo) are coupled 
coherently by a light of wave number kn . The coupling at the 
j- th site between the state 7 and 7 + 1 is denoted by The 
hyperhne states serve as an additional spatial dimension - the 
synthetic dimension. Since the phases of Q 7 depend on j, a 
particle hopping around a plaquette of this synthetic lattice 
picks up a phase equivalent to having a magnetic flux in the 
plaquette. 


controls their hopping from site to site. Further anal¬ 
ysis along these lines reveals, inter alia , a) the gauge 
held which induces a flavor-orbit coupling mitigates the 
“baryon breaking” effects of the Zeeman held, and b) 
the gauge held also induces a non local interaction, i.e., 
interaction between particles at different j sites. One 
crucial outcome is that under favourable circumstances, 
the SU{M ) singlet baryon (fly = 0), which is an ob¬ 
ject localized at a site j but extended along the synthetic 
dimension, is transformed into an M-body bound state 
that is extended in real space (along j ) which we dub 
as the “squished baryon”. This is demonstrated by ana¬ 
lytical arguments supported by detailed few body exact 
diagonalization calculations. This work not only suggests 
novel many-body phases (eg, “squished baryon” conden¬ 
sate) of these systems, but also suggests new opportuni¬ 
ties with SU{M) symmetric systems. For example, our 
mapping of the system to a set of particles with SU (M) 
gauge held brings out the possibility of using the M- 
component cold atom systems to simulate the SU{M) 
random flux model. m 


Methodology: Denoting the operators that creates a 
fermion[42 at site j with a hyperhne flavor 7 as Cj , 


the Hamiltonian is H = H t + Hq + Hu , with 


H t 

= ~t E (^+iTi.7 + h.c.) 

(1) 


j,7=1,M 


H n 

= E (^TF+x^+h. c.) 

(2) 


j, r y=l,M~ 1 


Hu 

= -j£ c L c L- c ,^ 

(3) 


J?77 / 


where t is the intersite hopping, U is the strength of the 
attractive SU(M) symmetric interaction. The couplings 
Qi, = £ 2 7 e -2/c ^ , where depend on details (see [37j 
[38]) of the system. 

A mapping gains further insights into the physics of 
the problem. Towards this end we introduce the notation 
where C j is the column (Cyi, Cj : 2 ,..., Cj : M) T - We in¬ 
troduce a local unitary transformation C j = VI j hj where 
VIj = Diag{e _z/c ? x -?}, 7 = 1,M, with kj = (7 — 1 )feg, and 
hj = (b jA ,...,b. ],m) T another set of fermionic oper¬ 
ators. Immediately, Hq = J2j LjOtbj where (D is a site 
independent Hermitian matrix 


0 = 


/ 0 nj 0 

f^l 0 f^2 

0 0 2 0 


0 \ 
0 


: : : ’• 1 

\o 0 ... n M -i 0 / 


(4) 


Diagonalization gives 0 = S<loS^ where <lo = Diagjco^}, 
£ = 1 is the diagonal matrix with eigenvalues 

and S is a unitary matrix. This results in Hq = 
oJojOGj, where S^tb j = ®j = {aj^} T , £ = 1,..., M is a 
new set of fermion operators. Clearly, Q = where 
= WjS is a unitary matrix. We now have 

H = -*52 (d+i u ]+i (L t'®i+ h - c -) + XbF®? + H u 

3 j 

(5) 

where Hu is operator defined in eqn. © written in terms 
of owing to its SU (. M ) invariance. We immediately 
see that in terms of the transformed states <Uj , the Hamil¬ 
tonian can be interpreted as that of particles in a flavour 
(() dependent potential (which is a SU{M ) Zeeman 
field), and whose hopping is influenced by a non-Abelian 

gauge field uJ +1 Uj = S^#S, (# = Diag je 2/e ^j) that 
produces flavor mixing up on hopping from site to site. 
The Zeeman field is determined solely by fi 7 , while the 
gauge field has a crucial additional dependence on the 
flux p/q. We have thus mapped the problem of synthetic 
dimensions to a system of M component fermions expe¬ 
riencing SU{M ) Zeeman and gauge fields and SU{M ) 
symmetric interactions. 

Induced Interactions: We now discuss the key out¬ 
come that arises in the analysis of eqn. Consider 
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FIG. 2. (Color online) Non-local interaction: Top panel 
shows the state of two fermions when t — 0 with M — 2 which 
has — —O and uj ^=2 — .$2. For Q < Q c , Q c = the 
usual baryon which is a bound state with both the f at a site 
being occupied (top left). When £2 > £2 C , baryon breaking 
effect of the Zeeman field sets in, and the lowest energy state 
is each of the two particles occupying £ = 1 state at arbitrary 
but distinct sites. At £2 = £2 C , the two states are degenerate. 
Arrows in the left bottom panel show the hopping pattern 
when t ^ 0 in the presence of a |-flux (- « §)• ^ the 

two particles of the broken baryon are in tlie neighbouring 
sites as shown, then this baryon can effectively hop on a dual 
lattice shown by crosses (bottom right) by hybridizing with 
the degenerate baryon (vertical shaded bond), gaining kinetic 
energy. This produces a net attractive interaction between 
particles at neighbouring site with £ = 1 and leads to “baryon 
squishing”, i. e., increasing non-local character of the bound 
state. 


the M = 2 system with p/q = 1/2, i.e., the so called \- 
flux case. The main idea is uncovered by looking at the 
rather unnatural limit of vanishing hopping t 0. The 
Zeeman field in this case is do = Diag{—£2,12}. When 
£2 <C E/, the ground state of the system is given by a 
M — 2 baryon with two particles localized at the same 
site (Fig. [2j top-left). The “baryon breaking” effect of 
the Zeeman field occurs when £2 exceeds 12 c = y (see 
Fig. [2j top-right). In the broken baryon state, the par¬ 
ticles both have £ = 1 and can be located at any two 
sites with distinct j. Turning on £, t > 0 with a |-flux 
now produces a hopping pattern of individual particles 
as shown by arrows in fig. [5] (bottom left). The hopping 
does not conserve the hyperfine flavour - this flavour- 
orbital coupling is a result of the gauge potential that 
depends on the flux. It is clear that the degeneracy of 
the broken baryon states is lifted by the hopping with 
flavor-orbit coupling - two particles with £ = 1 the on 
neighboring sites can gain energy by hybridizing with the 
degenerate baryon state (bound along the synthetic di¬ 
mension). The effect just discussed, therefore, induces 
a non local attractive interaction between particles with 
£ = 1 located on two neighbouring sites. Then net out¬ 
come of this is a “squished baryon” state that generically 


has a bound state character along the synthetic and real 
dimension. For the particular case just discussed, the 
state is a bound state of two particles that “resonates” 
between the vertical and horizontal bonds (as indicated 
bottom left figure of fig.[ 2 |, hopping on the “dual lattice” 
indicated by crosses in fig. [ 2 ] (bottom right). Indeed, as 
£2 >> 12 c , the bound state is primarily made of particles 
with ( = 1 - “fully” squished baryon which is a result of 
the attractive interaction between near-neighbour ( = 1 
states which is proportional to 2 q~u • 

Similar physics applies to a generic M. An important 
point to be noted is that the scale £ 2 C and the resulting 
“broken baryon” state depends on the details of the cou¬ 
pling f2 7 . For a given M and £2 7 , our arguments suggests 
that there are some special fluxes p/q that could more ef¬ 
fectively produce the effect of the non-local binding and 
baryon squishing just discussed. 

Exact diagonalization: We have numerically inves¬ 
tigated the few body (N p particle) physics in this sys¬ 
tem using N q sites of the physical lattice with periodic 
boundary conditions, resulting in a Hilbert space size of 
NqM Cn p • Our diagonalization scheme uses the transla¬ 
tional symmetry, with Q, the total center of mass mo¬ 
mentum, being the associated good quantum number. 
If overall ground state is found to occur with Q = Q g , 
then the binding energy is defined as E & = E g (Q g , U = 
0) — E g (Q g , E7), where E g (Q g , U = 0) is the ground state 
energy of the same system but without interactions be¬ 
tween the particles. We also study the properties of 
ground states by computing the moment of inertia along 
the x-direction, and average value for the synthetic coor¬ 
dinate £, 

^ = 1 vV<£( Aa: ^) 2 >’ ( 6 ) 

* 1>*2 

<c> = t/<Ec;>, (7) 

iv p i 

where is run over the particle labels and A Xi lr i 2 = 
( Xi t — X{ 2 ). An indicator of an TV^-particle bound state 
is a positive binding energy E^ > 0 and, I xx both of 
which are insensitive to N q (system size).|44] The quan¬ 
tity ((} provides a measure of squishing. For example, 
with N p = M, (C) = indicates the usual SU(M)- 

singlet baryon, while squishing can be deduced from a 
value of (£) < ( M +F . 

Results: While we choose the simplest case S2 7 = £2 
to illustrate the physical ideas, our calculations can be 
readily adapted for specific experimental systems. Fig. [3] 
shows the results for M = 2. In the absence of a flux 
p/q -A 0, the critical Zeeman field to break the baryon 
is Q c = \ (VU 2 + 1 6t 2 — 4 1) . The phase diagram in the 
p/q— £2 plane shown in fig.[3[a) and (b), show that this in¬ 
deed happens at p/q = 0. For larger £2, there is no bound 
state at p/q = 0. Notice, however, that for p/q = 1/2 
(-flux), the situation is entirely different. One sees that 
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FIG. 3. (Color online) M — 2 : (a) and (b) “Phase diagram” of two particles showing dependence (a) I XX: and (b) (£) on flux 
p/q and Q/U for t/U = 0.1. Insets show the type of bound state stabilized, (c) to (i) |-flux Eb , I xx and (£) are respectively 
shown in (c), (d) and (e). (f) shows the dependence of Eb on Q/U for t/U = 0.1, while (g) shows I xx and (0 for the same case. 
Dependence of Eb (h) and I xx , (£) on t/U at Q = Q c = U/2 . The dashed lines in (h) are results of analytical considerations at 
small and large jj. 
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FIG. 4. (Color online) M — 3 : (a) and (b) Phase diagram 3-particle phase diagram showing dependence (a) I xx , and (b) 
(C) on flux p/q and Q/U for t/U = 0.1 obtained with N q = 18. Insets show the type of bound state stabilized, (c) and (d), 
f-flux : Panels (c) and (d) show the dependence of the binding energy, and I xx and (0 on Q/U with t/U = 0.1. (e) to (h), 
p/q = 1/3: (e) and (f) show same quantities as (c) and (d) for the 1/3-flux case, (g) and (h) show the variation of Eb and 
Ixx, (C) as a function of t/U for Q = Q c = for the 1/3-flux case. 
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FIG. 5. (Color online) (Color online) Mm. 4 : (a) and (b) Phase diagram 4-particle phase diagram showing dependence 
(a) I xx , and (b) (() on flux p/q and Q/U for t/U = 0.1 obtained with N q = 8. Insets show the type of bound state stabilized, 
(c) and (d), |-flux: Panels (c) and (d) show the dependence of the binding energy, and I xx and (Q on Q/U with t/U = 0.1 
(e) and (f), p/q = 1/4: (e) and (f) show same quantities as (c) and (d) for the 1/4-flux case. 


the I xx remains finite with the increase in $2, and in fact 
(0 goes to unity. As shown in the insets of the phase 
diagram the baryon evolves to the squished baryon. We 
now investigate the \ flux case which has this interesting 
physics in greater depth. Figs. [3jc),(d) and (e) clearly 
demonstrate that for the —flux a bound state always ex¬ 
ists (except when t = 0) irrespective of a large Zeeman 
like field. This is a vivid example of the gauge field mit¬ 


igating the baryon breaking effects of the Zeeman field. 
Fig. [3jf) and (g) clearly demonstrate the squishing of 
the baryon by the Zeeman field, aided by the gauge field. 
Finally, fig. §h) and (i) are to discuss the case where 
= Q c . As shown from the analytic considerations (see 
fig-i, the binding energy of the squished baryon when 
t <C U is « 2 t as it involves the hybridization process 
discussed in fig. [2] Indeed, this result also corroborates 
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quantitatively with the earlier arguments. For example, 
the binding energy at small t is indeed found to be 21 . 

We now begin the discussion of the M m 3 case, whose 
results are shown in fig. [4j with a discussion of D c . When 
t = 0, D c = ^ with a peculiar feature. Three distinct 
states are degenerate at this value of D. There are the 
usual M = 3 barvon 127: . a completely broken baryon 
with three particles at different sites (“ 1 + 1 + 1 ”), and par¬ 
tially broken “ 2 + 1 ” baryon which has two particles at a 
given site with £ = 1 and 2 and the third particle at a dif¬ 
ferent site with £ = 1. Fig. |4ja) and (b) show the phase 
diagram for this case in the p/q-Q/U plane. Again the 
squishing effect of the Zeeman field aided by the non- 
Abelian gauge field is clearly seen. A crucial point is 
illustrated by the figs. |4jc)-(f); the (c) and (d) panels 
of the figure are for the case with a |-flux (t/U = 0.1 
), which show the squishing of the baryon continuously 
(most rapidly near ft c ) with increase of £1. However, the 
process does not go on forever, and at a value of Q some¬ 
what larger than D c , the baryon completely breaks up 
and there is no bound state. Here, therefore, the gauge 
field produced by the |-flux is unable to prevent the pair 
breaking effect. Most interestingly, the situation changes 
completely if one introduces a 1 /3-flux. As shown in 
fig -0 e) and (f), the squishing occurs smoothly, and in 
fact, we believe, that there is a bound state for all Q (we 
cannot verify this as binding energy becomes small with 
a concomitant large I xx ). Another interesting result is 
that at fi c , the binding energy for small t can be inferred 
to be proportional to t. This is due to the hybridization 
between the “ 2 + 1 ” baryon hybridizing with a “ 1 + 1 + 1 ” 
aided by the 1/3-gauge field (flavor-orbit coupling). In¬ 
deed, as shown in fig. Qg), this is in excellent agreement 
with the exact result. 

Results for M = 4 are shown in fig. [5j The novel as¬ 
pect here is the presence of two critical Zeeman fields 
fi c i and £l C 2 - When t = 0, the usual 4-baryon is desta¬ 
bilized to state with two 2 -baryons (each of which can 
be located at any site) at Q c i = ^=. At Cl C 2 = U this 
state is again broken into a l + l + l + l state where 
each particle can be at any site distinct from others with 
C = 1. With a |-flux, one sees a smooth transformation 
form the usual 4-baryon to a 2+2 baryon (bound state of 
2-baryons) - a vivid example of squishing. However, the 
|-flux is not able to mitigate the effects of the Zeeman 
field. In fact, near D c 2 the squished 2+2 baryon is bro¬ 
ken up. Remarkably, for a flux of 1/4, this transition is 
prevented, and our calculations suggest bound state for 
any D (checking this requires very large computational 
resources). At D c i, we can also show that the binding 
energy is proportional to t 2 (in-order to hybridize the 4- 
baryon and the 2+2 baryon); our exact calculations (not 
shown here) have borne this out. 

We now discuss the general criteria that are required 
to produce squishing (rather than breaking). Clearly, the 


non-abelian gauge field induced by the flux must be able 
to hybridize the two (or more) degenerate states that 
occur at the critical Zeeman fields. For example, the 
flavour-orbit coupling with a |-fhix does not hybridize 
the 2+2 state with the l + l + l + l state for M = 4, 
and hence the baryon is broken rather than squished. 
Knowing the details of fi 7 , one can choose appropriate 
flux to achieve this for a given system. 

In the many body setting, there is clearly a rich col¬ 
lection of phases and crossovers that can be explored. 
In particular, quasi-condensates of squished baryons are 
likely to hold interesting physics. The few body results 
developed here can be used as a guide for such studies 
particularly in the dilute limit. 

We conclude the paper by pointing out an interesting 
possibility to use this system to create a class Hamiltoni¬ 
ans call “random flux” models |4l . The idea is to intro¬ 
duce some randomness in D/, which in turn will make the 
gauge fields (eqn. (§) also random. If the lattice chosen 
is a square lattice, the Hamiltonian of the type (eqn. < 0 » 
realized will be similar to a “random flux” model. 
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